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THE THERMOSPHERIC TIDES 
AS AN 

OSCILLATOR CIRCUIT SYSTEM 
II. Volland 


ABSTRACT 

The in situ generation of neutral air waves within the thermosphere by an 
external heat source is calculated using an approximate approach where the 
density amplitude represents the electric current and the heat source replaces 
the electric voltage. The energy equation can be transformed in such a way that 
it formally is equivalent to the equation of an electric circuit system in series 
by neglecting the transmission quality of the waves. The " inductance” of the 
circuit is related to internal heat of the wave and to vertical heat convection, 
the "capacitance" is related to horizontal and vertical heat convection and the 
"resistance" is related to ion drag and vertical heat conduction. Numerical cal- 
culations depending on frequency, collision number and horizontal wave number 
agree sufficiently well with exact full wave calculations. They are applied to the 
diur*ial and semidiurnal tidal density variations, to the geomagnetic activity 
effect, the 27 day variation and the semiannual effect within thermospheric 
heights . 
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THE THERMOSPHERIC TIDES AS AN 
OSCILLATOR CIRCUIT SYSTEM 

X. INTRODUCTION 

Within the thermosphere the diurnal component of the tidal wave predominates 
the density variations. This tidal wave is generated by solar EUV heat input 
within the thermosphere on the one hand and by a diurnal tidal gravity wave from 
the lower atmosphere on the other hand (Volland, 1968). Above 300 km altitude 
the EUV generated tidal wave dominates. The observations show maximum 
density between 14° 0 and 15° 0 local time. The maximum is delayed by 2 to 3 
hours with respect to the solar EUV heat input which has its maximum at 
12 00 local time (see e.g., Priester et al. t 1967). It has been shown that this time 
delay is the natural response time of the thermospheric circuit system with re- 
spect to the generating heat source (Volland, 1968). 

We will give in this paper a simplified derivation of the diurnal tidal density 
wave within the thermosphere. We shall compare it with an electric network 
driven by an alternative voltage generator. We shall find complete equivalence 
between both systems. We therefore can define an atmospheric "condensor", 

1 ’inductor” and "resistor" equivalent to heat convection, internal energy and heat 
conduction and ion drag, respectively. The amplitude and phase relationship be- 
tween atmospheric "current" (density) and atmospheric "voltage" (heat source) 
and its dependence on the atmospheric "circuit elements" can be determined. 

This equivalence not only is valid for the diurnal period but generally for atmos- 
pheric disturbances of sinusodial form at thermospheric heights. 
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2. THE ELECTRIC CIRCUIT SYSTEM 


An electric circuit system in series consisting of inductance L, capacitance 
C and resistance R driven by an alternative voltage generator of voltage u (Bee 
Figure 1) has the well known relationship between electric current i and electric 
voltage u of 


i = Y (X) 

where u 0 is the amplitude and co is the angular frequency of the driving alterna- 
tive voltage generator , t is the time , 



is the impedance and 


( 2 ) 


$ 



(3) 


is the phase delay between current i and voltage u (see e.g., Glasford, 1965). 


Current and voltage are in phase if the frequency oj equals the eigenfrequency 
Q of the system: 



(4) 


The current i is delayed by 90° if R = 0 and co > Cl. i shows a phase advance of 
90° with respect to u if R = 0 and co < 0 . 
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3. THE ATMOSPHERIC CIRCUIT SYSTEM 


We now shall derive a completely equivalent formula for the diurnal tidal 
thermospheric density wave Ap generated by a heat input Aq. We consider the 
thermosphere as a circuit system in which energy is provided via direct EUV 
heat input Aq. Strictly speaking the thermosphere behaves like a transmission 
system of four port quality. Waves generated outside the thermosphere (the 
tidal wave from the lower atmosphere) propagate through the thermosphere with 
finite phase velocity and with finite attenuation rate, Waves generated at a cer- 
tain place within the thermosphere (e.g., by EUV heating) propagate out of the 
thermosphere and contribute to the diurnal density variation not only at the place 
of their generation but throughout the whole thermosphere. In the height range 
between 300 and 400 km, which we shall consider exclusively, the tidal wave 
from below still contributes significantly to the density variation. In our approxi- 
mation we neglect the transmission line quality of the thermosphere and consider 
only the in situ generation of the waves without taking into account their propa- 
gation characteristics. 

We consider a two dimensional model of an isothermal thermosphere and a 
heat input harmonic in local time r and exponentially decreasing with height z: 

Aq a Aq 0 e jc,J ' r “ (2 “ 2 o> /H (5) 

with 

cot - cot + k v . 

y J 
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o) is the angular frequency of the earth’s rotation, y is the coordinate along 
the equator (positive in east direction), z is the altitude, k y 58 1/a is the hori- 
zontal wave number (a is the earth’s radius), H is a scale height and z 0 is a 
reference height. From the equations of conservation of mass, momentum and 
energy and from the equation of! state we can derive the following expressions 
(see e.g., Volland, 1969a) 


BAw Po 

iwAp + p 0 + jk yl 0 o Av - -jj- Aw = 0 


(6a) 


jco/jpAv + v P q Av + jk Ap = 0 


(6b) 


3Ap 

TT + 8 = 0 


(6c) 


3Aw 3 / 3AT \ 

j wc v ^ 0 AT + jk yPo Av + p 0 («o- 37 -) = Aq 


(6d) 


Ap A p. AT 

mmmm » wmmmmrn 4* 

Po Po TT 


(6e) 


It is 

p 0 , T 0 , p 0 mean values of density, temperature and pressure 
A p r AT, Ap wave parameters of density, temperature and pressure 
Av, Aw horizontal wind (positive in east direction) and vertical wind 
v collision number between the ions and one neutral 
g magnitude of the gravitational acceleration force 
c v specific heat at constant volume 
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k q coefficient of mean heat conductivity 
H 0 scale height of the isothermal thermosphere. 

In the equation of conservation of vertical momentum (6c) we neglected the 
vertical velocity and in the equation of conservation of energy (Gd) we neglected 
the horizontal heal flux, because both terms are small compared with the vertical 
pressure gradient and the vertical heat flux, respectively. Since the wave propa- 
gation is considered to take place at low latitudes from west to east, the Coriolis 
force is small and can be neglected, too, in the equation of conservation of hori- 
zontal momentum (6b). Furthermore, in Equation (6b) we ignored molecular 
viscosity which gives rise to a term on the left side in Equation (6b) of 


- 7 ) 



( 7 ) 


where r) is the coefficient of molecular viscosity. 

From our approach given below we can estimate the magnitude of this term 
as compared with the ion drag term vyo 0 Av. At 350 km altitude using reasonable 
numerical data we find that the viscosity term is only a fraction of the ion drag 
term and can be considered to be included within an effective collision number y. 

The energy equation (6d) is in fact the first law of thermodynamics. The 

respective terms on the left hand side of Equation (6d) represent therefore the 

* 

change of the internal energy of the gas, horizontal and vertical heat convection 
(the work done by the gas) and vertical heat conduction considered as an ex- 
ternal heat sink, while the right hand side gives the change of the external 
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heat input. We shall make use of these expressions throughout the following 
sections. 


From an exact treatment of the complete system Equation (6) it can be shown 

BAT 

(Volland, 1968) that Ap, Ap and * 0 ar e nearly in phase and that their height 
dependence in the height range between 300 and 400 km can fairly well be de- 
scribed by exponential expressions: 


Ap 



o 

<3 

Ap 



^Pq | 

* 0 Ar 

> 


<3 

o 
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e - Hi 


( 8 ) 


Hj is a scale height and z 0 is the same reference height as in Equation (5). 

A p 0 , Ap 0 and k 0 ATq are related to this height. (For abb-v /i at ion we replaced 
the height derivative by a prime: 


B 




The absence of an imaginary term in the exponents of Equation (8) indicates the 
quasi -evanescent behavior of the tidal wave at those heights which is due to the 
dominant influence of heat conduction. 

■ *t 

i 

: 1 

Because we assume an isothermal thermosphere it is 

* 
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and 



( 10 ) 


The approximate nature of our approach becomes evident by the fact that the 
condition of an isothermal thermosphere (T 0 « const) r,nd of the exponential de- 
crease with altitude of the heat flux of the wave [k 0 AT* « exp (-z/faj) in Equa- 
tion (8)] can not be exactly fulfilled simultaneously in Equations (6). Moreover, 
within a realistic atmosphere scale height H 0 and collision number v are height 
dependent. Therefore, the scale height Hj of the wave parameters expected to 
be equal to the scale height H of the heat source has to be considered as an 
effective value slightly depending on altitude and on the parameters of the ther- 
mosphere. From an exact calculation one finds 

H 0 < Hj < H. 


Using Equations (8), (9) and (10) we can derive the following relationship 
between the different wave parameters from Equations (§)i 


Ap _ A p 
Po Pq 



, IN AP 


Av _ cxS A p 

yz 77 


c 


( 11 ) 



p o y 


c - velocity of sound 

y = ratio between the specific heats at constant pressure and at constant 
volume. 

Since the thermosphere is an open system, the external heat absorbed 
within a sheet at the altitude z 0 generates neutral air waves which re-radiate 
into the upper and into the lower half space, separated by z 0 . Therefore, within 
the upper half space z 2. z Q only half of the external energy Aq 0 is available to 
create the observed upgoing density wave. Thus, a factor 1/2 must be added to 
the external heat source Aq in Equation (6d) in order to meet such radiation 
condition. 


S 



Eliminating now the density amplitude A p from Equations (11) and 
this value into Equation (6d) yields 

Ze'^A p 0 = Aq 0 

with 

Z = Zei* = R + j («L - 


introducinj 


(12) 


(13) 


and 


R = R x + R 2 


a 2 k 2 c 4 v k q AT' 

i + 2 

7 2 (a> 2 + v 2 ) 


L 


2 (a - 1) c 2 
(7-1) 


C 


y 2 

2a 2 k y 2 c 4 



+ 



(14) 


Comparing Equations (12) and (1) we notice a complete equivalence, if we 
relate the density amplitude A p Q at the height z 0 to an electric current i and the 
heat input Aq 0 to an electric voltage u, and we can immediately define an atmos- 
pheric "inductor" L, "condensor" C and "resistor" R by Equations (14). 

From the Equation of energy (6d) we identify the "inductance" L as connected 
with the internal energy of the wave and with vertical heat convection, the "ca- 
pacitance" C as connected with horizontal and vertical heat convection and the 
"resistance" R as connected with icn drag and vertical heat conduction. The 
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be 


" capacitance” C increases with collision number v, with decreasing horizontal 
wave number k y and with decreasing frequency a>. The ion drag term Yi i of the 
"resistance” R first increases with v, reaches a maximum at v = co and then 
approaches zero for v 00 , it decreases with w and increases with k . The 
"inductance” L and the heat conduction term R 2 of R in our approximation re- 
main independent of co 9 v or k y . 

The eigenfrequency of the system according to Equation (4) is 


n = 


ak y c / 

2y y (a- 


(7 - 1) 

1) (1 + v 2 /a> 2 ) 


(15) 


However for k y f 0 Equation (15) because of the frequency dependence of 0 
is not identical with the effective eigenfrequency of the system for which the 
imaginary part of Z vanishes. This effective eigenfrequency is 


O) 


off 



2 kj 2 c J (y _ 

4y s (a - 1) 


(16) 


4. NUMERICAL RESULTS 

We shall give in this section numerical results of Equation (13). We shall 
compare these results with exact calculations in order to test the validity of 
Equation (13) and to find out some basic data about the relationship between 
density amplitude and heat input depending on frequency, collision number and 
horizontal wave number. We confine ourselves to a height of z 0 = 350 km and 
use as atmospheric parameters the data of the Jacchia-model (Jacchia, 1964) at 
an exospheric temperature of = 1000° K. 
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4a. Diurnal Density Variation and its Dependence on Collision Number v 

We first estimate the influence of the collision number v on the amplitude 
and phase relationship between diurnal EUV heat input and density. We use the 
data 

to = 7.28 x 10“ 5 sec’ 1 ; k y = 1. 5 x 10“ 4 km" 1 

7 = 1, 5 ; c = 0,84 km/sec 


H q ~ 53 km ; 

Hj = 74 km 

a = 1.4 ; 

and vary v from 10” 5 to 10“ 2 sec" 1 . 

’’impedance" Z and phase angle O [see Equation (13)] are plotted versus v as 
full lines. The dashed lines in Figure 2 have been calculated using the exact 
system of the dynamic equations, talcing into account the transmission charac- 
teristics of the thermosphere between 100 and 400 km altitude and the height 
dependence of v used in (Vollana, 1968). 

The number a = Hj/Hq in Equation (17) has been selected such that at z 0 = 
350 km and for v = 4 x 10" 4 sec" 1 the results of the exact calculations are in 
good agreement with the approximate relations in Equation (11). The number 
k q ATq / A p Q has been chosen in order to find optimal agreement within the whole 
range of collision number v and of frequency to. Actually both numbers slightly 


H - 100 km 


(17) 


«o AT o 
~ A ?o 


= 1.9 X 10’ 3 km 3 /sec 3 


The result is shown in Figure 2 where 
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depend on v and w. In our approximation we take them as constants * Therefore, 
we can not expect complete agreement between the exact calculations and our 
approximation. Nevertheless, as Figure 2 and the following figures indicate, the 
essential features of the density variation are well reproduced by our approxi- 
mation. 

Figure 2a shows that the "impedance” Z only slightly decreases with the 
collision number v . Therefore , the density amplitude generated by the EUV 
heat source does not depend very much on v . According to Jacchia (1964) the 
diurnal density amplitude at 350 km altitude and at moderate solar activity 
(F = 125) is 

| A/0,,1 = 3.3 x 10" 15 g/an 3 (18) 

Therefore the EUV heat input necessary to generate this wave amplitude at 
350 km altitude is 

| Aq 0 1 < Z|Ap 0 | “• 5 x 10“ 9 erg/cm 3 sec (19) 

The "smaller than" sign in Equation (19) takes account of the fact that the 
density amplitude of Equation (18) is in part generated by a tidal wave from 
below (Volland, 1968). 

The phase delay O of the density in Figure 2b varies from -75° to 70° with 
increasing collision number equivalent to a time advance and a time delay, re- 
spectively, with respect to the solar EUV heat source varying from -5 hours to 
4,7 hours. The observed time lag of about 2.5 hours in the density maximum 
therefore corresponds to a collision number of 
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V 


= 4 x 10"* 4 sec" 1 


( 20 ) 


which is in good agreement with the value derived from Dalgarno’s (1964) formula 
of v at the height of the maximum of the equatorial F2 layer. 


Figure 3a shows the v -dependence of the "circuit” elements coL, l/o> C, R x 
and R 2 [see Equation (14)] . In Figure 3b the real part 

R B Ri * R 2 (21) 


and the imaginary part 


of the complex impedance 


I = Cl)L 


coC 


( 22 ) 


Z = R + j I 

are plotted versus v. We notice that with increasing v , the influence of heat 
convection (the term X/coC) decreases, is compensated by the internal heat (the 
term oj L) at ci> = a> e{ { and finally disappears. The phase 0 in turn shifts from 



at 


over 



$ = 0 
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at 


/ 


a 2 k 2 c 2 (7 - 1) 


4 7 2 (a- 1) 


- co 2 


to 



at 

v -> 00 (23) 

Resistance R t due to ion drag is of the same order of magnitude than heat 
conduction term R 2 in the range v = 0 , Outside this range it decreases. 

4b. Semidiurnal Tides and Higher Harmonics 

The observations of the tidal density variation within the thermosphere 
clearly gives evidence of the predominance of the diurnal component (Jacchia 
and Slowey, 1966). This is in contrast to the situation within the lower atmos- 
phere, where the semidiurnal component predominates (see e.g., Siebert, 1961). 
The reason for the smallness of the diurnal component in the lower atmosphere 
in spite of its largest excitement by solar radiation as the driving force is due 
to the fact that this wave mode is of an evanescent type there (Siebert, 1961; 
Kato, 1966; Lindzen, 1967) so that its propagation is suppressed. 

Within the thermosphere the diurnal component changes into a propagation 
mode (Volland and Mayr, 1968), and its "impedance” Z becomes smaller than 



the respective values of the higher harmonics. In order to show this we repeated 
our calculations taking into account higher harmonics of the local time 


mr * eo m t + k m y (24) 

with 

- m co ; m - 1, 2, 3, 4 

k m - m k 

and plotted ’’impedance" Z and phase angle versus m in Figure 4. Since the 
diurnal component of the EUV heat source (exactly speaking the P u -component) 
is by far the largest component, and, moreover, Z in Figure 4a increases with 
m , the higher harmonics of the tidal wave are expected to be of no significance 
within the thermosphere, which is confirmed by the observations. 

4c. Geomagnetic Activity Effect, 27 Day Variation and Semiannual 
Variation of the Thermospheric Density 

The geomagnetic activity effect, the 27 day variation and the semiannual 
effect of the thermospheric density are nearly independent of solar local time 
(see e.g., Priester et al., 1967). Therefore, they can be treated by a one dimen- 
sional model (Volland, 1969b). In our approximation this implies the neglect of 
horizontal convection (k = 0). Then according to Equation (14) the influence of 
horizontal heat convection (1/wC) and of ion drag (R t ) disappears. Internal 
energy (the term co L, which includes a contribution of vertical heat convection) 
and vertical heat conduction (the term R 2 ) remain the only energy forms into 
which external solar EUV heat can be transformed. 
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In order to show the influence of frequency « on the excitement of a density 
wave we calculated Z and <I> as functions of eo using the same parameters as in 
Equation (17) (see Figure 5, full lines). The phase shifts from zero to about 60° 
within the frequency range considered and tends to become 00° for very high 
frequencies. We notice a slow increase of Z with oj, At periods greater than 
10 days the "impedance" Z remains nearly constant while the phase angle $ has 
dropped to zero. This follows from the decreasing influence of internal heat. It 
means that in this range of frequency the thermosphere behaves quasistatic ally 
and can be treated by the static equation of heat conduction. 

The dashed lines in Figure 5 are exact calculations of thermospheric dy- 
namics (Volland, 1969b). The four respective dashed lines in Figure 5 have been 
calculated taking into account a lower boundary of the heat input at 

z . = 100, 120, 150 and 200 km 

min 

altitude, respectively. Here we observe the interesting phenomen that the heat 
input within the lower thermosphere below 200 km altitude becomes of significant 
influence on the resulting values of Z and 0 at 350 km altitude for periods 
greater than one solar day. Our approximation is best in agreement with the 
150 or 200 km boundary data which means that the density wave, generated by 
the heat input between 100 and 150 km, is suppressed from propagating upwards 
at periods 

a 277* . 

At 0 - — £ 1 day . 

However, at higher periods the density wave can penetrate into the upper 
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thermosphere with increasing significance. Such phenomen can be called a 
thermospheric skin effect. 

The vertical bars in Figure 5b indicate the range of observed delay times 
of the geomagnetic activity effect with respect to the geomagnetic activity factor 
K [5 hours according to Roemer (1967a), 7 hours according to Jacchia et al. 

r 

(1966)1 , the delay time of the 27 day variation with respect to the excessive 
solar flux [1 day according to Roemer (1967b), 2.3 days according to MacDonald 
(1963)1 and of the semiannual effect with respect to the equinox [10 to 20 days 
according to Jacchia et al. (1968)] . In the case of the geomagnetic activity effect 
we assumed a quasiperiod of At 0 ~ 1.2 days which is an average impulse time 
of this effect. 

From a comparison between observations and calculations we conclude that 
during the 27 day variation the penetration of the solar EUV heat into the ther- 
mosphere stops at about 120 km altitude to be efficient. However, in the case of 
the semiannual effect most of the heat is deposited within 100 and 120 km altitude. 
This heat input is explained as arising from the dissipation of wave energy of 
the tidal wave from the lower atmosphere, the amplitude of which varies semi- 
annually due to changing prevailing winds within the lower atmosphere (Volland, 
1969b). Finally, the heat input which causes the geomagnetic activity effect also 
is deposited predominantly within the lower thermosphere. 

In this connection we notice from Figure 5 that the natural response time of 
the whole thermosphere with respect to diurnally varying heating is of the order 
of 5 hours, which is just the value found by Harris and Priester (1962) in their 
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one -dimensional model of the thermosphere without "second heat source”. 
Comparing I?igures 5 and 2 we see that horizontal heat correction is an essential 
part in the diurnal variation of the thermospheric density. 

5. CONCLUSION 

Using an approximate approach by neglecting the transmission quality of 
thermospheric waves the equations of the dynamics of the neutral thermosphere 
can be transformed in such a way that the equation of conservation of energy 
formally equals the equation of an electric circuit system. The "electric cur- 
rent" In this system is represented by the density amplitude, while the "electric 
voltage" is equivalent to the external heat input. The "inductance" of the system 
is related to internal heat and to vertical heat convection of the thermospheric 
wave, the "capacitance" is related to horizontal and vertical heat convection and 
the "resistance" is related to ion drag and vertical heat convection. 

In a numerical calculation depending on frequency, collision number and 
horizontal wave number the behavior of the different "circuit elements" is dis- 
cussed# The numerical results are compared with results from an exact full 
wave theory and show sufficient agreement for periods equal or smaller than one 
solar day. However for larger periods, waves generated within the lower ther- 
mosphere below 200 km can propagate upwards and significantly participate in 
the observed density variations. This behavior can be considered as a thermo- 
spheric skin effect. 

These calculations are applied to the diurnal and semidiurnal tidal density 

t 

variations within the thermosphere, to the geomagnetic activity affect, the 27 day 


viWatlon and the semiannual effect at thermospheric heights. The ratio between 
external heat source and density amplitude (the "impedance" of the circuit) 
and the phase delay between heat source and density maximum are determined 
and compared with observations. 
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FIGURE CAPTIONS 


Figure 1. 
Figure 2. 


Figure 3. 
Figure 4. 


Figure 5. 


Block diagram of an electric circuit system in series. 

Impedance Z (ratio between external heat source and density ampli- 
tude) and phase delay <I> of the density amplitude with respect to the 
external heat source of the diurnal thermospheric tidal density wave 
versus collision number v . Full lines have been calculated using 
the approximate approach, dashed lines are determined from an 
exact full wave calculation. 

fa) Circuit elements of the atmospheric circuit system. 

1 

(b) Real part R = R t + R 2 and imaginary part I = ojL - ^ of 
the complex impedance Z of the diurnal thermospheric tidal 
density wave versus collision number. 

Impedance Z and phase delay $ of the thermospheric tidal density 
wave depending on the domain number m . (m = 1 : fundamental or 
diurnal component ; m = 2 : semidiurnal component ; m > 2 : 
higher harmonics). Full lines have been calculated using the ap- 
proximate approach, dashed lines are determined from an exact full 
wave calculation. 

Impedance Z and phase delay $ of a one dimensional thermospheric 
model (k y = 0) versus angular frequency co (full lines). The upper 
scale gives the period (in days) of the waves. The dashed lines have 
been determined by an exact full wave calculation taking into account 



different lower boundaries z min of the external heat input, The ver- 
tical bars indicate ranges of observed phase delays of the density 
amplitude with respect to the supposed heat sources of the geomagnetic 
activity effect, the 27 day variation and of the semiannual effect. 



i (Ap 0 ) ► 

Figure 1, Block Diagram of on Electric Circuit System in Series 
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Figure 2. Impedance Z (ratio between external heat source and density amplitude) and phase delay $ of the density amplitude with 
respect to the external heat source of the diurnai thermospheric tidal density wave versus collision number v. Full lines have been 
calculated using the approximate approach, dashed lines are determined from an exact full wave calculation. 
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Figure 3. (a) Circuit elements of the atmospheric circuit system, (b) Real part R - R 1 + R 2 and imaginary part I 
of the complex impedance Z of the diurnal thermospheric tidal density wave versus collision number. 
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Figure 4. Impedance Z and phase delay $ of the thermospheric tidal density wave depending on 
the domain number m . (m = 1 : fundamental or diurnal component ; m = 2 i semidiurnal component; 
m > 2 : higher harmonics). Full lines have been calculated using the approximate approach, dashed 
lines are determined from an exact full wave calculation. 
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